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ABSTRACT 

In this paper we present a dust emission map of the starless core TMC-IC taken at 2100 /im. 
Along with maps at 160, 450, 850 and 1200 /xm, we study the dust emissivity spectral index from 
the (sub)millimeter spectral energy distribution, and find that it is close to the typically assumed 
value oi (3 — 2. We also map the dust temperature and column density in TMC-IC, and find that at 
the position of the dust peak [Ay ^ 50), the line-of-sight-averaged temperature is ~7 K. Employing 
simple Monte Carlo modeling, we show that the data are consistent with a constant value for the 
emissivity spectral index over the whole map of TMC-IC. 

Subject headings: stars: formation; ISM: dust, extinction 
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1. INTRODUCTION 

Starless cores are often identified through maps of 
their dust emission at (sub)millimeter wavelengths, and 
recent surveys of nearby molecular clouds have in- 
creased the number of such ob jects co nsiderably (e.g. 
HatchcU et al. 2007; En och et al. ' 2008; K auffmann et al.l 
2008; .Simpson et al..,2008f) . Analvsis of the properties of 
starless cores is often restricted by sparse sampling of 
their spectral energy distribution (SED), uncertain dust 
emissivities and assumptions about their geometries. As 
a result, the measured (sub)millimcter fluxes are often 
used in conjunction with assumed values for the dust 
temperature and/or emissivity to determine core masses, 
though the systematic errors can be substantial (a factor 
of a few). With sufficient spectral coverage and knowl- 
ege of source geometry, it should be possible to map the 
column density, temperature and dust emissivity spec- 
tral index in a starless core, though this has yet to be 
accomplished. 

One of the best-studied starless cores is TMC-IC, 
located in the Taurus m olecular cl oud at an approxi- 
mate distance of 140 pc ({Torres e~ al. 2009). Previous 
studies of the dust emission from TMC-IC at 450, 850 
and 1200 /im have determined that it is cold {Td ~6 
K) and dense {n ^10^ cm~^) at its center, and be- 
comes progressively l e ss dense an d warmer at larger radii 
(jSchnee et ail l200'7al: ISchnee fc Goodman 2005.1 . The 
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dust-derived mass of TMC-IC is a factor of a few larger 
than its virial mass, and self-abs orbed N2H+ spectra 
show evidence for sub-sonic infall (|Schnee et al.ll2007bl : 
ISchnee fc Goodman|[2005l ) . The chemistry of TMC-IC is 
characterized by depletion of CO and its i sotopologues, 
as well as depletion of CS and N2H+ (jSchnee et al.l 
[2007b). 

Previously reported mass and temperature profiles for 
TMC-IC are estimated assuming a constant value for th e 
emissivity spectral index (/3 — 1.8) (jSchnee et al.ll2007al ). 
However, a wide range of values for (3 in other regions 
have been measured or calculated on theoretical grounds. 
For instance, in the diffuse interstellar medium (ISM), 
an emissivity spectral i ndex of /3 = 2 fits the typical far- 
infrar ed spectrum well (|Draine et al.ll2007tlDraine fc Lei 
[T98l . while grain g rowth in protostellar disk s results 
in values of /3 ~ 1 (jBeckwith fc Sargentlll99l( l. Since 
starless cores are in an intermediate stage between the 
diffuse ISM and protostars, one might expect that the 
dust in a starless core would have an emissivity spectral 
index in the range 1 < /? < 2, although the timescale 
of grain growth at densities nn ~ 10^ cm^'^ is several 
hund red times the f ree-fall time of a typical starless core 
(Cha krabarti fc Mc Kcc 2005). On the other hand, the 
growth of icy mantles on dust grains could steepen the 
slope of the dust SED, and this has been cited as the 
reason for high values of the emissivit y spectral index 
in parts of the Orion Ridge (3 2.5) (iLis et all [19981) 
and in Sgr B2(N) (/3 ~ 3.7) (IKuan et al.lll996D . In the 
cold and dense evironment in TMC-IC, the dust grains 
could form similar icy mantles, which might drive the 
emissivity spectral index to values /3 > 2. 

In addition to changes in /3, grain growth by either 
coagulation or accretion of an icy mantle will change 
the opacity of the dust grains, k. For instance, ev- 
idence that dust opacities are higher in dense cores 
than in the diffuse ISM can be found by comparing gas 
temp erature and dust emission profiles (jKeto fc Casellil 
[200l or by comp aring dust and gas emission profiles 
(iKeto et al l I2004D . From models of dust coagulation 



for conditions that are pla usible for a starless core, 
lOssenkopf fc Hennind (|1994f ) show that the dust opac- 
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ity at 1.3 mm can be a factor of five higher than in 
a diffuse environmen t. Using 160 yum data for Taurus, 
ITerebev et al.l (|2009[ ) measure an opacity of kibo fim — 
0.23 ± 0.046 cm^g~^ in the cold cloud (Ay ranging from 
0.4 to 4.0). The op acity is 2.6 times greater tha n the 
diffuse ISM opacity (jWeingartner h Draind [20031 ) . and 
provides e vidence for higher opacity v alues in dense re- 
gions (c.f. lOssenkopf fc H ennine) 1 1994( 1. 

Laboratory studies of possible dust grain materials 
have measured the dust emission and absorption prop- 
erties at far-infrared an d longer wavelengths. For in- 
stance, lAannestadl (|1975l ) measured the absorption of sil- 
icate grains and found that /3 ~ 2 for grains without an 
ice mantle, and /3 ~ 3 for grains with an ice mantle. 
The steepest spectral indices calculated by lAannestadl 
(fl975.) were for fused quartz and olivine with ice man- 
tles, both of which had / 3 = 3.5, clos e to th e extremely 
steep /3 = 3.7 reported bv lKuan et al.l ()1996l ). In a study 
of the millimeter-wave absorptio n of silicate gr a ins in 
the temperature range 1.2 - 30 K, lAgladze et all (IT996h 
found that the emissivity spectral index depends on the 
dust temperature and for some materials can vary be- 
tween 1.5 < /3 < 2.5, with the largest value for /3 at a 
temperature ~ 10 K. A temperature- dependent emis- 
sivity spectral index was also found by iMennella et al.l 
(|1998l ) for cosmic dust analog grains in the tempera- 
ture range 24 < < 295 K and a wavelength range 
20 [iva. < A < 2 mm, with a steeper SED at lower tem- 
peratures. The range of emissivity spectral indices is 
1.2 < (3 < 2.4 for the various grain types considered by 
IMennella et al.l ()1998[ ). In a study of amorphous silicates 
in the temperature range 10 < Td < 300 K and wave- 
length range 0.1 < A < 2 mm. lBoudet et aTl ([2005) also 
report an anticorrelation between Td and /3, with values 
of the emissivity spectral index between 1.5 < /3 < 2.5. 
Based on these laboratory measurements, one might ex- 
pect that in starless cores li ke TMC-IC, wh ich have tem- 
peratures close to 10 K (Sc hnee et al.ll2009l ). the emissiv- 
ity spectral index would be close to /? ~ 2.5. 

Observations of nearby star-forming regions suggest 
an anticorrelation between Td and /3 similar to that 
measured in the laboratory. For instance, observa- 
tions of the Orion, p Ophiuchi and Taurus molecular 
clouds at far-infrared and (sub)millimeter wavelengths 
by Dupac et al. (2003) show that the emissivity spectral 
index may vary from 2.5 > /3 > 1.0 as the dust temper- 
ature rises from 10 < T^; < 80 K. Similarly, in a study 
of the far-infrared and (sub)millimeter S EDs of luminous 
infrared galaxies. I Yang &: PhillipsI (|2007| ) report that the 
emissivity spectral index varies from 2.4 > /3 > 0.9 as 
the dust temperature varies from 30 < Td < 60 K. In 
a study of the starless core L1498, S hirlev et all ([2005h 
find an emissivity spectral index of /3 = 2.44 ±0.62 for an 
assumed dust temperature ofTd = 10.5±0.5K, consis- 
tent with the hypothesis that the cold dust in a starless 
core has a steep spectral index. 

However, variations in the dust properties along the 
line-of-sight and noise in the obser ved maps can mimic an 
anticorrelation between Td and /3 (jShettv et al.ll2009al lbi) . 
Furthermore, at high densities and low temperatures, 
grain coagulation will decrease P while the possible T^- 
P relation will increase the emissivity spectral index. 
Therefore, it is not at all clear that one should expect the 



emissivity spectral index to vary significantly in a star- 
less core. Observations of dust emission and extinction 
in the dark cloud B68 s how no evidence for variations in 
the dust properties (Bi anchi et al.l[2003f l. suggesting that 
to some extent the various pressures on /? may cancel out. 

In this paper we present maps of the starless core 
TMC-IC at 160, 450, 850, 1200 and 2100 pm. By spa- 
tially resolving the emission from TMC-IC over a wide 
range of wavelengths, we then map the column density, 
temperature and emissivity spectral index of the dust, 
which has typically been difhcult to do with previous 
datasets. 

2. OBSERVATIONS 

The 160, 450, 850 and 1200 /xm maps have been pre- 
sente d or described in previous p apers (ITerebev et alJ 
lOOg"; 'Schnee & Goodman' 'Schnee et aL| |2007a|; 

Di Francesco ct al. 2008; Kauffmann ct al. 200^, so here 
we briefly review the observations and data reduction, 
and refer the reader to the earlier papers for more details. 
The 2100 /zm map has not been previously presented, so 
it is discussed in detail below. 

2.1. Spitzer 

The Spitzer data near T MC-IC derive from the Tau- 
rus Spitzer Legacy Survey (jPadgett et al.ll2009l ). a large 
44 square degree map of the Taurus star-forming region 
at seven wave l ength s ranging from 3.6 to 160 fim. See 
iPadgett et all ([2009) for an overview of the survey pa- 
rameters and data processing. The MIPS data utilize 
the fast scan mode with 3 seconds integration time per 
frame and 5 visits per pixel. The MIPS 160 /xm data ob- 
tained in fast scan mode mapping does not have enough 
redundancy to fill in completely all the data gaps due to 
a dead readout and other effects. To mitigate these data 
gaps and to preserve the diffuse emission as much as pos- 
sible, a 3 pixel by 3 pixel (15" per pixel) median filter is 
applied to the image. The median filter introduces some 
image smoothing, increasing the effective angular resolu- 
tion from 40" to 1'. In some cases, filtering reduces the 
surface brightness by 10-15%, though the reduction in 
the map of TMC-IC is only 2%. 

Analysis of the 160 /im data in Taurus shows there are 
two main cloud components: emission from the extended 
cold cloud (14.2 ± 0.4 K , Ay < 3), and emission from 
compact cold co res (Td ~10 K, high Ay) (jTerebev et al.l 
120091 : IfTagev et al. 2009). Since our present focus is the 
high extinction regions sampled by (chopped) millime- 
ter data, here we further process the 160 /im data to 
remove the spatially extended cloud emi ssion. Follow- 
ing the procedure in ITerebev et al.l (|2009[ ) we use IRAS 
100 fim images to model and subtract the 14.2 K cloud 
emission. The resulting 160 /im image traces the high 
extinction ridge and looks very similar to the millimeter 
continuum emission. The median flux removed from orig- 
inal 160 /im map is 96 MJy/sr, or about 90% of the total 
flux. Uncertainties in the temperature of the extended 
cloud component results in a 10% uncertainty in the final 
160 fim map. Given t he 12% absolute flux un certainty of 
MIPS 160 fim maps (jStansberrv et al.ll2007f l. we conser- 
vatively estimate the absolute uncertainty in the 160 fim 
"cold core" map to be 20%. 

2.2. SCUBA 
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Although in previous papers on TMC-IC we used 
sub-millimeter maps observed and reduced by us 
(jSchnee fc Goodman! [200l iSchnee et al.ll2007aD . in this 
paper we use maps of TMC-IC at 45 and 850 jum from 
the new SCUBA Legacy Catalogues ()Di Francesco et all 
[20081 . We choose to use the new maps be cause they in- 
clude all of the data taken with SCUBA (|Holland et al.l 
|1999[) at all epochs and use the most current calibra- 
tion and reduction methods. The noise in the maps are 
132 and 21 mJy beam~^ at 450 and 850 fim, respec- 
tively, and have resolutions of appoximately 9" and 14". 
The uncertainties in the absolute flux calibration (pri- 
marily due to beam-shape uncertainty and fluctuations 
of opacity above the telescope during observations and 
calibration) are 20% at 850 jum and 5 0% at 450 /im, as 
described in iDi Francesco et all (|2008f ) . We set the zero 
point of the Legacy maps by determining the median flux 
density in two regions off the core but nearby and sub- 
tracting these values (-0.32 and -0.035 Jy beam^^) from 
the 450 and 850 /im maps. 

2.3. MAMBO 

We observed the 1200 /xm continuum emission from 
TMC-IC in the autumns of 2002 a nd 2003 using the 
MAMBO-2 array (jKrevsa et al.lll999f ) on the IRAM 30 m 
telescope on Pico Veleta (Spain), producing a map with 
10.7" resolution. Data were taken with on-the-fly map- 
ping, with the telescope sub-reflector chopping in az- 
imuth by 60"-70", and the resultant im age was recon- 
structed using the lEmerson et al.l (|1979[ ) algorithm iter- 
atively to properly reproduce large-scale emission. The 
noise in map is 3 mJy beam~^, and the flux calibration 
uncertainty, as derived from the rms of calibrator obser- 
vations across pooled observing sessions and the uncer- 
tainty in the intrinsic calibrator fluxes, is about 10%. For 
further deta ils on the 1200 /xm ob serva tions and data re- 
duction, see ISchnee et al.l ()2007aD and iKauffmann et all 
l|2008, ). 

2.4. Bolocam 

We observed the 2100 /im continuum emission from 
TMC-IC in the autumn of 2007 using Bolocam on the 
Caltech Submillimeter Observatory (CSO). Data were 
collected in raster scan mode at a scan speed of 240"s~^, 
with a step size between rasters of 11". Each raster was 
20' in length, and each of the 105 total scans were sep- 
arated by 11". The resultant map has a coverage of 
~20'x20'with 60" resolution. Three sets of observations 
contain rasters parallel to RA, and the other three con- 
tain rasters parallel to Dec. Data reduction and flux cal- 
ibration, accomplished through observatio ns of Uranus 
and N eptune, was done as described in I Savers et al.l 
(l2009l ). 

Observations at approximately 2 mm are dominated by 
sky noise. Although the Bolocam pipeline PCA cleaning 
algorithm efficiently removes sky noise, it also removes 
some of the source flux. This decreases the peak in- 
tensity, introduces negative bowls around bright sources, 
and reduces the brightness on scales larger than 1' for ex- 
tended sources. We recover this lost flux using a mapping 
algorithm that iteratively subtracts a source model from 
the real da ta (for more d escription and performance eval- 
uation, see lEnoch et al.li2006,) . This technique is similar 



to CLEAN (IHogbomlll97"l ISchwarzlll97l . but works in 
the image plane, and the real data are used for the source 
model rather than, e.g., a Gaussian component model. 
The first source model is formed from the original map, 
including only regions with signal-to-noise greater than 
1.5. That source model is subtracted from the original 
timestream data, which is then re-cleaned and mapped 
without the source flux. This process is iterated until the 
source model stabilizes, 10 iterations in this case. TMC- 
IC is a faint source, so the correction is relatively small, 
with the source flux density increasing by up to 30% af- 
ter iterative mapping. The noise is the 2100 /im map is 
8 mJy beam"^, and the uncertainty in the absolute flux 
calibration is ^20%. 

3. ANALYSIS 

Here we estimate the temperature, column density and 
emissivity spectral index in TMC-IC using our flux den- 
sity maps at 160, 450, 850, 1200 and 2100 /tm. Under 
the assumptions that the dust is isothermal and char- 
acterized by a single emissivity spectral index, the flux 
density per beam in each map can be described as a mod- 
ified blackbody, given by, 

= ^}B^{Td)^i^^J.mHNH2 (1) 

where 

Biy{Td) - — (2) 

eyi\){hv/ kid) — 1 

and 

"''^"^^"(230^)' 

In equation ([T|), S^, is the flux density per beam, Vl 
is the solid angle of the beam, Bv(Td) is the blackbody 
emission from the dust at temperature T^, /i = 2.8 is the 
mean molecular weight of interstellar material in a molec- 
ular cloud per hydrogen molecule, mH is the mass of the 
hydrogen atom, is the column density of hydrogen 
molecules, and a gas-to-dust ratio of 100 is assumed. In 
equation = 0.009 cm^ g~^ is the emissivity of 

the dust grains at a gas density of 10^ cm~^ covered 
by a thin ice mantle at 230 GHz (jOssenkopf fc Hennind 
119941 Column 6 of Table 1) and 13 is the emissivity spec- 
tral index of the dust. Although Equations ([T]-[ni) are 
often used to described the observed emission from star- 
less cores, such observations average over a range of dust 
grain properties. Furthermore, the emissivity (^230) is 
uncertain by roughly a factor of 2. If we had chosen 
a different value for K230, this would affect the derived 
column density, but would have no impact on the de- 
rived temperature and emissivity spectral index. We do 
not consider variations in K230 in this paper, but instead 
save this analysis for a future analysis. 

3.1. Color Correction 

The fluxes given in the dust emission maps are de- 
rived from the total energy detected convolved with the 
wavelength-dependent instrumental response and an as- 
sumed SED of the emitting source. The assumed SEDs 
in the Spitzer, SCUBA, MAMBO and Bolocam maps are 
that of a 10,000 K blackbody, a flat-spectrum source, a 
flat-spectrum source and a blackbody in the Rayleigh- 
Jeans regime, respectively. Because we are observing 
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emission from dust, with an assumed SED of a modified 
blackbody, the quoted flux at the nominal wavelength of 
each map ought to be corrected. However, the relatively 
narrow bandpasses of the detectors results in a modest 
correction factor. 

For each map, we use the instrumental response and 
our best-guess dust SED of a modified blackbody with 
Td = 10 K and /3 = 2 to find the "reference" wavelength 
of the observations. The integral of the convolution of 
the instrumental response and the dust SED are equal 
above and below this reference wavelength. 

/ SxR{X)dX= / SxRiX)dX (4) 

where S\ is the flux emitted by the dust at wavelength A 
and R{X) is the instrumental response at wavelength A. 

The reference wavelength, used in our calculations of 
the dust emission properties, are close to the nominal 
wavelength of each map, with the biggest difference in 
the MAMBO map because of its wide bandpass. In all 
maps, the flux uncertainties are larger than the correc- 
tion from using the reference wavelength in our calcula- 
tions rather than the nominal wavelength. The derived 
reference wavelengths are given in Table [1] 

3.2. Single SED 

One way to estimate the dust temperature, column 
density and emissivity spectral index is to fit modified 
blackbody SEDs (given by Equations [1] - [3]) to the five 
fluxes (at 160, 450, 850, 1200 and 2100 yum) at each po- 
sition. Although the resolutions of the maps at 450, 850 
and 1200 fim are fairly similar (from 7.5" - 14"), the 60" 
resolution of the 160 and 2100 /im maps is not a good 
match to the others. Here we handle the mismatched 
beam sizes by making an average SED from the five flux 
maps to derive the best- fit dust properties (Td, and 
/3). 

We first smooth the 450, 850 and 1200 /im maps to the 
resolution of the 160 and 2100 /im maps by convolution 
with a two-dimensional Gaussian beam. We then cal- 
culate the average flux in those pixels brighter than 3cr 
in the 160, 850, 1200, and 2100 /xm maps, and brighter 
than 2(7 in the 450 /im map (which is noisier than the 
other maps) . The area with significant emission in all five 
maps is shown in Figure [TJ The noise level in the maps is 
determined in regions off-source but near TMC-IC (see 
Figure [1]). We then used a minimization routine to 
fit the five measured fluxes with a modified blackbody 
SED, finding the values of T^, and /3. 

In order to determine how noise and calibration uncer- 
tainties affect the accuracy of the derived dust properties, 
we run a simple Monte Carlo simulation. In each trial, 
Gaussian random noise of the amount observed is added 
to each flux map, and then each map is scaled by a mul- 
tiplicative factor / = 1.0-1-5, where 5 is randomly chosen 
from a normal distribution of mean zero and standard 
deviation equal to the absolute calibration uncertainty 
in each map. Each map is modified by a different scale 
factor /, and this process was repeated 5000 times. We 
find that the best-fit dust properties for the average SED 
over TMC-IC are /3 = 2.3 ± 0.4, Td = 8.9 ± 1.0, and 
= 10.2x10^1 ± 2.1x1021. The observed fluxes and 



best-fit SED are plotted in Figure [5] and the fluxes and 
uncertainties are given in Table [T] 

A fundamental flaw in this analysis is that equations 
[1] - |3] are only accurate for isothermal dust character- 
ized by a single emissivity spectral index, but TMC- 
IC is known to be warme r at its edges than at the 
center (jSchnee et al.l l200"7ah . and it is certainly possi- 
ble that 3 also changes with in the starless core. As 
shown bv lShettv et al.l (|2009al ). the incorrect assumption 
of isothermal dust drives an anti-correlation between the 
derived dust temperature and emissivity spectral index. 
A spurious anti-correlation between Td an d may also 
arise d ue to noise in the emission maps (jShettv et al.l 
l2009airbl ). An anti-correlation between /? and Td is seen 
in the Monte Carlo simulation we used to derive the dust 
properties, as shown in Figure [2l In the following sec- 
tion we attempt to mitigate the effects of variations in 
the dust properties on our analysis. 

3.3. Maps 

The dust properties (Td, Nh^ and (3) in TMC-IC may 
vary both in the plane of the sky and along the line of 
sight. The line-of-sight variations in the dust properties 
cannot be determined without knowledge of the geome- 
try of the core, but variations in the plane of the sky are 
more tractable. Here we make maps of TMC-IC, first of 
the temperature and column density while holding the 
emissivity spectral index constant, and then of all three 
quantities. 

3.3.1. Maps with Constant (3 

We first regrid the 450, 850 and 1200 /im maps to the 
resolution of the 850 /im map, 14". We then determine 
the noise in the maps in the same regions selected in 
Section 13.21 and calculate the dust temperature and col- 
umn density in those pixels brighter than 3f7 in the 850 
and 1200 /im maps and brighter than 2a in the 450 /im 
map, using equations [T] - [3] for an assume d value of the 
emissiv ity spectral index, as was done in iSchnee et al.l 
(l2007a| ). Temperature and column density maps are 
made for each value of the emissivity spectral index in 
the range 1.0 < /3 < 4.0 in steps of 0.1. To account 
for the effects of noise and calibration uncertainties, we 
use the same Monte Carlo technique described in Section 
l3.2l to add Gaussian random noise to each pixel and ad- 
just each map by a multiplicative scale factor, and this 
process is repeated 100 times for each value of /3. 

To determine which value of (3 is most appropriate for 
TMC-IC, we then use the derived Td and Nh maps to 
predict the fluxes at 160, 450, 850, 1200 and 2100 /xm. 
The predicted 160 and 2100 /im fluxes are convolved with 
a 60" Gaussian beam to match the resolution of the ob- 
served maps. We find that the 160 and 2100 /im maps 
are best reproduced by an emissivity spectral index of 
p = 2.2 ± 0.6. The observed, predicted and residual flux 
maps for [3 = 2.2 are shown in Figure [3l and the de- 
rived temperature and column density maps are shown 
in Figure SI The column density is expressed in units of 
Ay, using t he conversion facto r Nii2 = 9.4 x lO^'' cm^^ 
(Ay /mag) (jBohhn et al.llT978l ). 

3.3.2. Maps with Variable l3 

In section 13.3.11 we assumed that f3 is constant across 
TMC-IC to make 14" maps of the dust temperature 
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and column density in TMC-IC. Here we make two- 
dimensional maps of the dust temperature, column den- 
sity and emissivity spectral index in TMC-IC, with the 
tradeoff that the resultant maps have coarser resolution 
(60"). 

We first regrid the 450, 850 and 1200 /im maps to the 
60" resolution of the 160 and 2100 /im maps. For each 
pixel brighter than 3ct, we then calculate Td, and (3 
using equations [l] - [S] To estimate the uncertainties in 
this method due to noise and absolute calibration errors, 
wc used the same Monte Carlo techniques as in sections 
13.21 and 13.3.11 with 1000 trials. In each trial there were 
about 15 pixels with sufficiently high signal to noise to 
calculate the dust properties, and within the maps there 
is a clear anti-correlation between Td and /3, as shown in 
Figure m The median emissivity spectral index and dust 
temperature derived from the high signal to noise points 
in all trials are /3 = 2.0 ± 0.6 and Td = 9.4 ± 1.9, similar 
to that derived from the composite SED of TMC-IC in 
section 13.21 

4. DISCUSSION 



In Section 13.3.21 we found that there is an anti- 
correlation between the dust temperature and emissivity 
spectral index in the two dimensional maps of TMC-IC. 
Although this pheno menon has been rep orted in other 
star-forming regions (Pup ae et al.l |2003f) and in other 



galaxies (Yang & Phillips 20071 ) . this is not a known fea- 
ture of starless cores. The assumption that the temper- 
ature along each line-of-sight is constant, even though 
the temperature profile in TMC-IC is very likely to be 
colder at higher densities, can lead to the mistaken con- 
clusion that there is an anti-correlation between Td and f3 
(jShettv et al.ll2009al lbl). Noise in the dust emission maps 
can also drive a spurious T^ — (3 anti-correlation. Here we 
investigate whether the dust temperature and emissivity 
spectral index are truely anti-correlated, and conclude 
that this is not necessarily the case. We then determine 
what constant value of (3 best characterizes TMC-IC. 

4.1. Testing the Td — (3 Anti- correlation 

To determine the effects of temperature variations 
along the line of sight and noise on the dust proper- 
ties derived from flux maps at 160, 450, 850, 1200 and 
2100 /im, we consider a spherical model of a starless core 
with plausible properties. The density in the core is given 

by 

= 1 + (r7.o)- 

where r is the radius, uq = 10^ cm~'^ and tq — 0.005 pc. 
The temperature in the core is given by 



Td(r) =ro(H-r/ri) 



(6) 



where Tg = 6.4 K and ri = 0.15 pc. The emissivity 
spectral index is constant in the core at a value (3 = 2.2. 
This model core is colder and denser at the center, as 
expected for an externally heated starless core. We do 
not claim that these are the temperature and density 
profiles for TMC-IC, but there is a general agreement 
between the shape of the model SED and that of TMC- 
IC. The density and temperature profiles are shown in 
Figure [5l as is the total emergent SED. 



For each line of sight passing through the model core, 
we determine the volume of intersection between a cylin- 
der and each spherical shell of the core. Given the density 
and temperature distributions from Equations [5] and [6l 
the flux from each layer is given by 



Si, 



B,y{Td)HynrVr 



(7) 



where Vr is the volume of intersection and d is the dis- 
tance to the core. The two-dimensional flux distribution 
is calculated by summing the flux from each layer along 
each line of sight. For the sake of convenience, the total 
flux density of the model core, summed over all radii, 
is scaled by a common factor such that the peak of the 
SED has a value of unity. 

We "observe" the model core by deriving two- 
dimensional flux maps at 160, 450, 850, 1200 and 
2100 /im and adding Gaussian random noise equal to 
5% of the peak flux in each map. Using equations [T] - 
[3l we then derive two-dimensional maps of the column 
density, dust temperature and emissivity spectral index, 
assuming that these quantities are constant along the line 
of sight (as was done in Section r3. 3. 2p . The derived col- 
umn density map peaks at the center of the core, where 
the temperature is seen to decrease, as one would expect 
from the input density and temperature profiles (see Fig- 
ure [5]). The derived emissivity spectral index map covers 
a wide range of values, mostly between 1.0 < (3 < 3.5, 
even though the input emissivity spectral index was held 
constant. The median derived emissivity spectral index 
is 2.0 ± 0.4, which is close to the input value of /3 = 2.2. 

Furthermore, the derived temperature and emissivity 
spectral index are seen to be anti-correlated. The com- 
bination of noise and temperature variations along the 
line of sight have been shown t o drive an a nti-correlation 
between these two quantities ([Shettv et a l. 2009a, b), so 
it is not surprising that this effect is seen here. Inter- 
estingly, the portion of the Td — (3 parameter space cov- 
ered by the model starless core overlaps significantly with 
that of TMC-IC (see Figures [S^il and p. From this 
we conclude that there is no definitive evidence that the 
emissivity spectral index in TMC-IC varies with tem- 
perature, and that the more reliable estimates of column 
density and temperature will come from fitting the ob- 
served fluxes while holding the emissivity spectral index 
constant. 

4.2. Average Emissivity Spectral Index in TMC-IC 

Estimates of the emissivity spectral index in TMC-IC 
derived from its average SED (see Section [3?2| and from 
two dimensional maps using the 450, 850 and 1200 /im 
maps to predict the lower resolution 160 and 2100 /im 
maps (see Section EXJ]), f3 = 2.3 ±0.4 and /3 = 2.2 ±0.6, 
are in good agreement with each other. Given that an 
analysis of the flux maps from an idealized core, with 
linc-of-sight temperature variations and noise, gives the 
correct average emissivity spectral index (within errors) , 
we suggest that the value of the emissivity spectral index 
in TMC-IC is well-represented by /3 ~ 2.2 ± 0.5. This 
is consistent with the emissivity spectra l index in the 
diffus e ISM of /3 = 2 (|Draine et all [20071 : IDraine fc Lei 
119841) and with the emissivity spectral index measured 
in the starless core L1498 (/3 = 2.44 ± 0.62. , Shirlev et all 
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120051) . However, due to the large uncertainty in (3 in this 
study, we cannot rule out the possibility that the emissiv- 
ity spectral index in TMC-IC could also be intermediate 
between the shallower SEDs seen in circumstellar disks 
and the SED of the ISM. 

When deriving the dust temperature and mass of 
prestellar cores from sparsely-sampled SEDs, it is com- 
mon to assume that the emissivity spect ral index of the 
dust is equal to 2 fe.g. iKirk et al.l 120071: Simpso n et aT] 
[200l . Observations of L1498 (jShirlev et al.ii2005[ ) and 
our observations of TMC-IC are consistent with /3 = 2, 
though the best-fit values of the emissivity spectral in- 
dex in both cores are slightly steeper, which would re- 
sult in lower temperatures and higher masses. However, 
with so few measurements of P in starless cores, there 
is not sufficient information with which to re-examine 
the dust properties in previously published studies. Up- 
coming surveys of the Gould Belt molecular clouds using 
Herschef and SCUBA-2^ promise to make statistically 
significant measurements of the emissivity spectral index 
in starless cores. 

5. SUMMARY 

In this paper we have analyzed the dust emission at 
160, 450, 850, 1200 and 2100 /im from TMC-IC. Our 
previous studies of this starless core have mapped the 
dust temperature and column density. Here, using new 
160 and 2100 /im maps and an improved reduction of 
the 450 and 850 /im maps, we are able to estimate the 
emissivity spectral index and make an improved analysis 
of the temperature and column density in TMC-IC. 

We find that the spectral index of the dust emission, 
as calculated by a variety of related methods, is in the 

^ http:/ /starformation-herschel. iap.fr/gouldbelt 
^ http:/ /www.jach.hawaii.odu/JCMT/surveys/gb/ 



range 1.7 < /? < 2.7 This is consistent with the value of 
the emissivity spectral index often assumed in studies of 
starless cores, and is also consistent with laboratory mea- 
surements and observations of the behavior of dust grains 
at low temperatures (T^ < 10 K). However, there is still 
a relatively large uncertainty in the emissivity spectral 
index. 

Using the new estimate of the emissivity spectral in- 
dex in TMC-IC {(3 =2.2), we construct temperature and 
column density maps, which reach values of ~ 7 if 
and Ay — 50 at the position of the dust peak, as shown 
in Figure m We find no conclusive evidence that (3 varies 
within the core, though line-of-sight temperature varia- 
tions and noise in the flux maps leave open the possibility 
that the emissivity spectral index is anti-correlated with 
dust temperature. 



We thank our referee, Yancy Shirley, for comments 
that improved the clarity of this paper. SS acknowl- 
edges support from the Owens Valley Radio Observatory, 
which is supported by the National Science Foundation 
through grant AST 05-40399. JS was partially supported 
by a NASA Postdoctoral Program Fellowship. Support 
was provided to ME by NASA through the Spitzer Space 
Telescope Fellowship Program. The JCMT is operated 
by the Joint Astronomy Centre on behalf of the Par- 
ticle Physics and Astronomy Research Council of the 
United Kingdom, the Netherlands Organisation for Sci- 
entific Research, and the National Research Council of 
Canada. IRAM is supported by INSU/CNRS (France), 
MPG (Germany), and IGN (Spain). The CSO is sup- 
ported by the NSF fund under contract AST 02-29008. 



REFERENCES 



Aannestad, P. A. 1975, ApJ, 200, 30 

Agladze, N. I., Sievers, A. J., Jones, S. A., Burlitch, J. M., & 

Beckwith, S. V. W. 1996, ApJ, 462, 1026 
Beckwith, S. V. W., & Sargent, A. I. 1991, ApJ, 381, 250 
Bianchi, S., Davies, J. I., Alton, P. B., Gerin, M., & Casoli, F. 

2000, A&A, 353, L13 
Bianchi, S., Gongalves, J., Albrecht, M., Caselli, P., Chini, R., 

Gain, D., & Walmsley, M. 2003, A&A, 399, L43 
Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 
Boudet, N., Mutschke, H., Nayral, C., Jager, C., Bernard, J. -P., 

Henning, T., & Meny, C. 2005, ApJ, 633, 272 
Chakrabarti, S., & McKee, C. F. 2005, ApJ, 631, 792 
Di Francesco, J., Johnstone, D., Kirk, H., MacKenzie, T., & 

Ledwosinska, E. 2008, ApJS, 175, 277 
Draine, B. T., et al. 2007, ApJ, 663, 866 
Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 
Dupac, X., et al. 2003, A&A, 404, Lll 

Emerson, D. T., Klein, U., & Haslam, C. G. T. 1979, A&A, 76, 92 
Enoch, M. L., Evans, N. J., II, Sargent, A. I., Glenn, J., 

Rosolowsky, E., & Myers, P. 2008, ApJ, 684, 1240 
Enoch, M. L., et al. 2006, ApJ, 638, 293 
Flagey, N., et al. 2009, arXiv:0906.5025 

Hatchell, J., Fuller, G. A., Richer, J. S., Harries, T. J., & Ladd, 

E. F. 2007, A&A, 468, 1009 
Hogbom, J. A. 1974, A&AS, 15, 417 
Holland, W. S., et al. 1999, MNRAS, 303, 659 
Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, N. J., II, & 

Lee, C. W. 2008, A&A, 487, 993 
Keto, E., Rybicki, G. B., Bergin, E. A., & Plume, R. 2004, ApJ, 

613, 355 



Keto, E., & Caselli, P. 2008, ApJ, 683, 238 

Kirk, J. M., Ward-Thompson, D., & Andre, P. 2007, MNRAS, 
375, 843 

Kreysa, E., et al. 1999, Infrared Physics and Technology, 40, 191 
Kuan, Y.-J., Mehringer, D. M., & Snyder, L. E. 1996, ApJ, 459, 
619 

Lis, D. C, Serabyn, E., Keene, J., Dowell, C. D., Benford, D. J., 
PhiUips, T. G., Hunter, T. R., & Wang, N. 1998, ApJ, 509, 299 

Mennella, V., Brucato, J. R., Colangeli, L., Palumbo, P., 
Rotundi, A., & Bussoletti, E. 1998, ApJ, 496, 1058 

Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 

Padgett, D. L. et al. 2009, in preparation 

Pierce-Price, D., et al. 2000, ApJ, 545, L121 

Sayers, J., et al. 2009, ApJ, 690, 1597 

Schnee, S., Rosolowsky, E., Foster, J., Enoch, M., & Sargent, A. 

2009, ApJ, 691, 1754 
Schnee, S., Caselli, P., Goodman, A., Arce, H. G., 

Ballesteros-Paredes, J., & Kuchibhotla, K. 2007, ApJ, 671, 1839 
Schnee, S., Kauffmann, J., Goodman, A., & Bertoldi, F. 2007, 

ApJ, 657, 838 
Schnee, S., & Goodman, A. 2005, ApJ, 624, 254 
Schwarz, U. J. 1978, A&A, 65, 345 

Shetty, R., Kauffmann, J., Schnee, S., Goodman, A. A., & 

Ercolano, B. 2009, ApJ, 696, 2234 
Shetty, R., Kauffmann, J., Schnee, S., & Goodman, A. A. 2009, 

ApJ, 696, 676 

Shirley, Y. L., Nordhaus, M. K., Grcevich, J. M., Evans, N. J., II, 
Rawlings, J. M. C, & Tatematsu, K. 2005, ApJ, 632, 982 

Simpson, R. J., Nutter, D., & Ward-Thompson, D. 2008, 
MNRAS, 391, 205 



Emissivity Spectral Index in TMC-IC 



7 




300 200 100 -100-200-300 
RA offset (orcsec) 



Fig. 1. — Flux maps of TMC-IC at 160, 450, 850, 1200 and 2100 fim, in units of mjy beam~^, smoothed with the 2100 /jm map's beam. 
The white contour shows the region that is above Sit in the 160, 850, 1200 and 2100 fira maps and above 2a in the 450 ^m map, and is the 
region used for determining the dust properties. The white boxes show the regions over which the noise in the maps was calculated. The 
beam size of each map is shown in the bottom left corners of each panel. 
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Fig. 2.— (left) The fluxes of TMC-IC at 160, 450, 850, 1200 and 2100 ^lm, in a 60" Gaussian beam, averaged over the high signal-to-noise 
portion of the maps, shown with la error bars. The dotted line shows the best-flt modified blac kbod y spectrum, which has column density 
N[j^ = 10.2x10^^, dust temperature T^j = 8.9 and emissivity spectral index /3 = 2.3. See section r3.2l for details. (rigr/ii) Filled contours show 
the density of — /3 values derived in the Monte Carlo simulation of the average TMC-IC fluxes, with the open blue contours containing 
50% and 95% of all points. The yellow circle shows the position of the median dust temperature and emissivity spectral index, used to 
calculate the SED in the left side of this figure. The fluxes and uncertainties are given in Table [T] 



TABLE 1 
TMC-IC Fluxes 



Nominal Wavelength 


Reference Wavelength 


Average Flux 


Flux Uncertainty 






(Jy/beam-'- ) 




160 


159 


2.0 


20% 


450 


456 


5.6 


50% 


850 


856 


1.1 


20% 


1200 


1090 


0.70 


10% 


2100 


2103 


0.050 


20% 



in a 60'' Gaussian beam 
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Fig. 3. — (left column) The observed maps of TMC-IC at 160, 450, 850, 1200 and 2100 /xm, shown in units of mjy beam"^. (center 
column) The predicted maps of TMC-IC, based on the derived temperature and column density maps (see Figure|4}, with /3 =2.2. (right 
column) The absolute value of the residuals. See section FS.S.ll for details. 
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Fig. 4. — [top row) Column density and dust t emper ature maps, expressed in units of Ay and K, derived from 14" maps of TMC-IC at 
450, 850 and 1200 fim using /3 =2.2. See section [3]3!T] for details, {center row) Column density, dus t tem perature and emissivity spectral 
index maps of TMC-IC, derived from 60" maps at 160, 450, 850, 1200 and 2100 fim. See section [TX2l for details, (bottom row) Filled 
contours show the density of T^j — /3 values, with the open blue contours containing 50% and 95% of all points. The yellow circle shows the 
position of the median dust temperature (T^ = 9.4) and emissivity spectral index (/3 = 2.0). The red crosses show the values of /3 and 
derived from each independent position in the 60" flux density maps with high signal to noise. 
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Fig. 5. — {top left and center) Density and temperature profiles of a spherical model of a starless core, (top right) The SED of the model 
starless core with /3 = 2.2, normalized such that the SED peak is at unity. The TMC-IC observations are also plotted for reference, and 
scaled by a common factor such that the 450 fim flux from TMC-IC is equal to that of the model core, {center row) Column density, 
temperature and emissivity spectral index maps that would be derived from the noisy maps of the model starless core, {bottom row) Filled 
contours show the density of T^j — /3 values, with the open blue contours containing 50% and 95% of all points. The yellow circle shows the 
position of the median dust temperature and emissivity spectral index. The anti-correlation between and /3 is similar to that observed 
in TMC-IC (see Figuregll. 



